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Abstract. In models, nucleation mode ( 1 nm < D p < 10 nm) 
particle microphysics can be represented explicitly with 
aerosol microphysical processes or can be parameterized 
to obtain the growth and survival of nuclei to the model’s 
lower size boundary. This study investigates how the repre- 
sentation of nucleation mode microphysics impacts aerosol 
number predictions in the TwO-Moment Aerosol Sectional 
(TOMAS) aerosol microphysics model running with the 
GISS GCM Il-prime by varying its lowest diameter bound- 
ary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm bound- 
ary simulates the nucleation mode particles with fully re- 
solved microphysical processes, while the model with the 
10 nm and 3 nm boundaries uses a nucleation mode dynam- 
ics parameterization to account for the growth of nucleated 
particles to 10 nm and 3 nm, respectively. We also investigate 
the impact of the time step for aerosol microphysical pro- 
cesses (a 10 min versus a 1 h time step) to aerosol number 
predictions in the TOMAS models with explicit dynamics 
for the nucleation mode particles (i.e., 3 nm and 1 nm bound- 
ary). The model with the explicit microphysics (i.e., 1 nm 
boundary) with the 10 min time step is used as a numerical 
benchmark simulation to estimate biases caused by varying 
the lower size cutoff and the time step. Different represen- 
tations of the nucleation mode have a significant effect on 
the formation rate of particles larger than 10 nm from nucle- 
ated particles (J io) and the burdens and lifetimes ofultrafine- 
mode (10 nm < Dp < 70 nm) particles but have less impact 
on the burdens and lifetimes of CCN-sized particles. The 
models using parameterized microphysics (i.e., 10 nm and 


3 nm boundaries) result in higher J io and shorter coagula- 
tion lifetimes of ultrafine-mode particles than the model with 
explicit dynamics (i.e., 1 nm boundary). The spatial distribu- 
tions of CN10 (D p > 10 nm) and CCN(0.2 %) (i.e., CCN con- 
centrations at 0.2 % supersaturation) are moderately affected, 
especially CN10 predictions above ~700hPa where nucle- 
ation contributes most strongly to CN10 concentrations. The 
lowermost-layer CN10 is substantially improved with the 
3 nm boundary (compared to 10 nm) in most areas. The over- 
prediction in CN10 with the 3 nm and 10 nm boundaries can 
be explained by the overprediction of J \ o or J 3 with the pa- 
rameterized microphysics, possibly due to the instantaneous 
growth rate assumption in the survival and growth param- 
eterization. The errors in CN10 predictions are sensitive to 
the choice of the lower size boundary but not to the choice 
of the time step applied to the microphysical processes. The 
spatial distribution of CCN(0.2 %) with the 3 nm boundary is 
almost identical to that with the 1 nm boundary, but that with 
the lOnm boundary can differ more than 10-40% in some 
areas. We found that the deviation in the 10 nm simulations 
is partly due to the longer time step (i.e., 1 h time step used 
in the 10 nm simulations compared to 10 min time step used 
in the benchmark simulations), but, even with the same time 
step, the lOnm cutoff showed noticeably higher errors than 
the 3 nm cutoff. In conclusion, we generally recommend us- 
ing a lower diameter boundary of 3 nm for studies focused 
on aerosol indirect effects but down to 1 nm boundary for 
studies focused on CN10 predictions or nucleation. 
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1 Introduction 

Atmospheric aerosols perturb Earth’s energy balance by scat- 
tering and absorbing solar radiation, known as the aerosol 
direct effect, and modifying cloud microphysical properties 
by acting as cloud condensation nuclei (CCN), known as the 
aerosol indirect effect. The impact of anthropogenic aerosols 
on global radiative fluxes is significant but highly uncertain 
due to the aerosol indirect effect, which is the most uncertain 
of anthropogenic climate forcings. One of the uncertainties in 
the estimates of the aerosol indirect effect stems from chal- 
lenges in predicting global distributions of CCN. At a fixed 
supersaturation, the ability of a particle to act as a CCN is de- 
termined by the particle size and chemical composition and is 
described relatively well by Kohler theory (Seinfeld and Pan- 
dis, 1998). The large uncertainty in CCN prediction is mostly 
from the estimates of particle number concentrations at CCN 
sizes. Atmospheric particles, including CCN-sized particles, 
can be emitted directly from sources (i.e., primary emission) 
or can be formed through the microphysical growth of nucle- 
ated particles (the formation of ~ 1 nm particles from con- 
densable vapors). Uncertainty in these source rates can con- 
tribute to large uncertainties in CCN prediction (Pierce and 
Adams, 2009c). 

Regarding nucleation, there are several widely suggested 
theories: binary (H 2 O-H 2 SO 4 ) nucleation, ternary (H 2 O- 
H 2 SO 4 -NH 3 ) nucleation, and ion-induced nucleation (Na- 
pari et ah, 2002; Vehkamaki et ah, 2002; Yu et ah, 2008). In 
these theories, H 2 SO 4 , NH 3 , H 2 O and molecular ions clus- 
ter together, and when they reach a cluster diameter around 
1 nm, they become stable and will not evaporate (Kulmala 
et ah, 2004). The rates at which these stable ~ 1 nm parti- 
cles are formed vary by orders of magnitude between differ- 
ent nucleation theories and parameterizations under the same 
conditions (Kulmala et ah, 2004). Uncertainty in the nucle- 
ation rates can contribute to large uncertainties in CCN pre- 
diction if the growth of nucleated particles to CCN sizes is 
significant. Several studies have investigated the sensitivity 
of CCN predictions to nucleation rates (Kristjansson et ah, 
2008; Makkonen et ah, 2009; Merikanto et ah, 2009; Pierce 
and Adams, 2009b; Wang and Penner, 2009; Yu and Luo, 
2009). The sensitivity of CCN particles to nucleation rates 
will depend on other atmospheric conditions, including pri- 
mary emission rates and the amount of condensable vapor. 
However, most studies show less than 20 % changes in CCN 
concentrations at 0.2 % supersaturation (hereafter defined as 
CCN(0.2%)) in the boundary layer, resulting from uncer- 
tainties in the boundary layer nucleation rate (e.g., binary 
vs. ternary nucleation schemes in Pierce and Adams (2009b) 
and with vs. without activation nucleation in the boundary 
layer with binary nucleation above the boundary layer in 
Merikanto et ah, 2009). Yu and Luo (2009) present more 
than a factor of two larger impact of nucleation to bound- 
ary layer CCN(0.4%) (CCN concentrations at 0.4% su- 
persaturation) because it includes the contribution from the 


nucleation above the boundary layer (compared to the bound- 
ary layer nucleation only in other studies). Their reported im- 
pact may be large because it labels sulfate particles nucleated 
in plumes as nucleation, whereas other studies count them as 
primary particles. Although these studies have examined the 
effect of nucleation on CCN concentrations, there has been 
less study of how to numerically simulate nucleation prop- 
erly in a global aerosol model. 

Various global models with online sectional aerosol micro- 
physics have chosen different cutoffs for the smallest sim- 
ulated aerosol size. Several models have used a lower cut- 
off of 1 nm (Yu and Luo, 2009; Snow-Kropla et ah, 2011). 
The choice of 1 nm is logical because the critical cluster 
size is ~ 1 nm; thus the entire size distribution is explic- 
itly simulated. Other models have used lower size cutoffs 
of 3 nm (Spracklen et ah, 2005, 2008) and 10 nm (Adams 
and Seinfeld, 2002; Trivitayanurak et ah, 2008; Lee et ah, 
2009; Pierce and Adams, 2009c, b). The rationale for these 
larger cutoffs was several-fold: ( 1 ) the smallest particles have 
the shortest lifetimes with respect to coagulation, and thus 
their explicit simulation adds a disproportionate computa- 
tional burden to the model; ( 2 ) they are only present sporad- 
ically in the atmosphere; and (3) until recently atmospheric 
measurements of aerosol number were generally for parti- 
cles above 3 or 10 nm depending on the instrument. To avoid 
the high computational burden required to simulate the nu- 
cleation mode particles explicitly, a nucleation mode param- 
eterization is commonly used in global aerosol models. Ker- 
minen et ah (2004a) proposed a parameterization to account 
for the growth of fresh nuclei to larger sizes (hereafter re- 
ferred to as the Kerminen parameterization). The Kerminen 
parameterization takes into account condensational growth 
of fresh nuclei and their coagulational scavenging by pre- 
existing particles. However, it does not account for coagu- 
lational growth of nucleated particles (i.e., self-coagulation 
among nucleated particles), which may be important under 
high-nucleation events. We should mention that the most up- 
dated nucleation mode dynamics parameterization includes 
the effect of the nuclei self-coagulation (Anttila et ah, 2010) 
and the elimination of the slightly inaccurate coagulation 
sink (Lehtinen et ah, 2007) that are missing in the current ver- 
sion of GISS-TOMAS. Nonetheless, the Kerminen parame- 
terization avoids simulating explicit nucleation mode dynam- 
ics and reduces the computational demand in a global model. 
To our knowledge, on the global scale, how the application 
of this parameterization at various lowest size boundaries im- 
pacts the aerosol number budgets has not been tested against 
an explicit aerosol microphysics model. 

In this paper, we investigate the impact of the lowest 
size boundary on aerosol budgets and evaluate the Ker- 
minen parameterization with different lower size bound- 
aries against explicit nucleation mode dynamics, using the 
TOMAS aerosol microphysics module in the Goddard In- 
stitute for Space Studies General Circulation Model (GISS 
GCM) Il-prime. With a size cutoff smaller than 10 nm, a time 
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step in the aerosol microphysics module shorter than 1 h may 
be preferred due to the short coagulational lifetime of the nu- 
cleation mode particles (Kerminen et ah, 2004b). Thus, we 
also investigate the impact of the choice of time step in the 
TOMAS algorithm on aerosol number budgets. Sections 2 
and 3 provide the description of the model and the setup for 
the simulations used for this work. Section 4 presents how 
the model aerosol number budgets are affected by the lowest 
size limit and a time step applied to the aerosol microphysical 
processes, and Sect. 5 is the conclusions. 


2 GISS-TOMAS description 

The TwO-Moment Aerosol Sectional (TOMAS) aerosol mi- 
crophysics model (Adams and Seinfeld, 2002) has been im- 
plemented into the climate model of GISS GCM II-prime, 
referred to as the “GISS-TOMAS” model (Lee and Adams, 
2010). The GISS GCM II-prime has horizontal grid dimen- 
sions of 4° latitude and 5° longitude, with nine vertical sigma 
layers including the stratosphere to the lOhPa level (Hansen 
et al., 1983); the pressure levels of the boundaries for the nine 
layers are 984, 934, 854, 720, 550, 390, 255, 150, 70, and 
10 hPa. Only 1-2 layers are in the stratosphere, so this model 
essentially treats only tropospheric aerosol. A detailed de- 
scription of the GISS GCM is found in Hansen et al. (1983). 
The parameterizations of convective and stratiform clouds 
are updated by Del Genio and Yao (1993) and Del Genio et 
al. (1996), respectively. Chemical tracers are advected every 
hour by the model winds using a quadratic upstream scheme 
(Prather, 1986); heat and moisture are advected with a simi- 
lar scheme. The time step for tracer processes in the GCM is 
1 h, but the TOMAS aerosol microphysical processes use an 
internal, adaptive time step. 

The TOMAS aerosol microphysics model uses a sectional 
approach that represents the aerosol size distribution by pre- 
dicting the amount of aerosol in several size sections or 
“bins”. TOMAS tracks two moments of the aerosol size dis- 
tribution in each size bin: total aerosol number and mass. To- 
tal mass is decomposed into several aerosol species, allow- 
ing prediction of the size-resolved aerosol composition. The 
original configuration of the model has 30 size sections (de- 
noted TOMAS-30), with the lower boundary of the small- 
est size bin being 10 -21 kg dry mass, and each successive 
boundary has twice the mass of the previous boundary. This 
provides a size distribution that ranges approximately from 
10 nm to 10 pm in dry diameter, depending on aerosol den- 
sity (see Fig. 1). TOMAS uses a moving sectional approach 
to treat water uptake; changes in water mass do not move 
particles between sections. Adams and Seinfeld (2002) pro- 
vide a detailed description of the TOMAS model. The model 
tracks nine quantities for each size bin: sulfate mass, sea salt 
mass, mass of externally mixed elemental carbon (EC), mass 
of internally mixed EC (mixed with all other species), mass 
of hydrophobic organic matter (OM), mass of hydrophilic 
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Fig. 1. Configurations of three TOMAS models used in this study 
and the number of simulations used for each model. Note that the 
TOMAS-36 and TOMAS-30 models use the nucleation mode dy- 
namics parameterization (i.e., Kerminen parameterization) to ac- 
count for the growth of nucleated particles to the model’s lowest 
size boundary. 


OM, mass of mineral dust, mass of ammonium, and mass of 
aerosol water and the number of aerosol particles in that size 
section. In addition, the model tracks one bulk aerosol-phase 
species, methanesulfonic acid (MSA), and six bulk gas-phase 
species: H2O2, SO2, dimethylsulfide (DMS), H2SO4, ammo- 
nia (NH3), and a lumped gas-phase tracer representing oxi- 
dized organic vapors that can form secondary organic aerosol 
(SOA). Gas-phase H2SO4 is assumed to be in pseudo-steady- 
state equilibrium between chemical production and conden- 
sational/nucleation losses (Pierce and Adams, 2009a). Water 
uptake by sulfate, sea salt, and hydrophilic OM is accounted 
for in the model. 

The TOMAS model used in this work includes previ- 
ously developed modules for sulfate (Adams and Seinfeld, 
2002), sea salt (Pierce and Adams, 2006), and carbonaceous 
aerosols (Pierce et al., 2007), and mineral dust (Lee et al., 
2009). Emissions of primary particles are treated as in pre- 
vious TOMAS studies and are briefly summarized here. The 
primary sulfate is emitted in two modes: fifteen percent of 
the mass is emitted to the first mode (number median diam- 
eter (NMD) of 10 nm; geometric standard deviation (GSD) 
of 1.6), and the rest is emitted to the second mode (NMD of 
70 nm; GSD of 2) (Adams and Seinfeld, 2003). The size dis- 
tribution of biofuel and biomass burning carbonaceous par- 
ticles is assumed to have a NMD of 100 nm and a GSD of 
2. Fossil fuel carbonaceous emission is assumed to have two 
modes based on Ban- Weiss et al. (2010). About 1.8 % of the 
mass (64 % of the total number) is emitted to a smaller mode 
(NMD of 17.5 nm and GSD of 1.6) and the rest is emitted to a 
larger mode (NMD of 60 nm and GSD of 1.9). Sea salt emis- 
sion parameterization is based on Clarke et al. (2006), and 
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mineral dust emission parameterization is based on Ginoux 
et al. (2001) and Marticorena and Bergametti (1995). The 
TOMAS model has been evaluated with ground-level mea- 
surement number and mass concentrations, deposition fluxes, 
and remote sensing observations (Adams and Seinfeld 2002; 
Pierce and Adams 2006; Pierce et al., 2007; Lee et al., 2009; 
Lee and Adams, 2010). In addition, the TOMAS coagulation 
and condensation algorithms have been evaluated against an- 
alytical solutions and have shown excellent agreement (Jung 
et al., 2006; Lee and Adams, 2012). 

The descriptions of dry deposition and wet deposition are 
available in Adams and Seinfeld (2002) and Lee et al. (2009). 
Briefly, dry deposition uses the series resistance approach 
that treats size-dependent gravitational settling of particles 
and a size-dependent resistance in the quasi-laminar sub- 
layer. Wet deposition occurs in large-scale (stratiform) and 
convective clouds and only in warm clouds. For in-cloud 
scavenging, the model uses modified Kohler theory to de- 
termine whether or not a particle is a CCN. The large-scale 
and convective clouds assume to have supersaturations of 
0.2% and 1.0%, respectively. The activation calculation is 
performed in each grid cell and time step with precipitation 
based on instantaneous aerosol size and composition. CCN 
concentrations presented in this paper are computed based 
on the monthly mean aerosol size and chemical composition 
fields (Pierce et al., 2007; Lee et al., 2009 for more details). 
The particles larger than the activation diameter at 0.2 % su- 
persaturation is counted as CCN(0.2 %). When the activation 
diameter falls in between size boundaries, we perform in- 
terpolation to determine what fraction of particles in a size 
section activate. 


3 Representations of nucleation mode microphysics 

As shown in Fig. 1, the TOMAS-30 model uses the Kerminen 
parameterization to account for the growth of nucleated parti- 
cles up to the first size bin, a diameter of 10 nm, instead of ex- 
plicitly simulating aerosol microphysics below 10 nm (Pierce 
and Adams, 2009b). The Kerminen parameterization pre- 
dicts the formation rate of particles at the model’s lower size 
boundary (i.e., 10 nm in the standard TOMAS-30 model) that 
occurs due to the nucleation rate, the condensational growth 
rate of the nucleated clusters by sulfuric acid and SOA pre- 
cursor gases, and the coagulational losses of fresh nuclei 
with existing particles larger than the desired size (Kerminen 
et al., 2004a). The version of the parameterization used in 
this work does not consider the growth of the nuclei through 
self-coagulation and thus may lead to an underprediction of 
the 10 nm particle (or other desired size) formation rate in 
the presence of high new-particle formation rates (Kermi- 
nen et al., 2004a). Furthermore, Pierce and Adams (2009c) 
explain that the Kerminen parameterization assumes a con- 
stant growth rate from the critical cluster size to the model’s 
lower size boundary based on the instantaneous growth rate 


when the particles were nucleated. If the true growth rate de- 
creases during the time it takes a critical cluster to grow to 
the model’s lower size boundary of lOnm, then the Kermi- 
nen parameterization would overpredict the J\o rate. 

To find out how the model-predicted global aerosol con- 
centrations depend on the lower size cutoff (size at which 
the Kenninen parameterization predicts the formation rate of 
particles), three model configurations are used in this paper 
(shown in Fig. 1). TOMAS-30 is our original model config- 
uration that has lOnm as its lowest size boundary. The Ker- 
minen parameterization is used to predict the formation rate 
of 10 nm particles from growth of nuclei. The TOMAS-36 
model extends the lowest size boundary down to 3 nm, and 
the Kerminen parameterization is applied to predict the for- 
mation rate of 3 nm particles. The TOMAS-40 model has the 
lowest size boundary extended to 1 nm. Explicit aerosol mi- 
crophysics (i.e., the TOMAS model itself) is used to simu- 
late the nucleation mode dynamics in this case except that 
the Kerminen parameterization is still applied to get the for- 
mation rate of 1 nm particles when the diameter of critical 
nuclei particles is smaller than 1 nm. 

When the explicit microphysics is used for the nucleation 
mode particles (i.e., TOMAS-36 and TOMAS-40 models), it 
may be necessary to use a shorter time step in the aerosol mi- 
crophysics module than 1 h as used in the TOMAS-30 model 
due to the fast microphysics of the nucleation mode particles 
(Kerminen et al., 2004b). To find out how the choice of a 
time step for aerosol microphysical processes affects aerosol 
number budgets, a 10 min time step is used for TOMAS-36 
and TOMAS-40 models in one set of simulations, and the 
1 h time step (i.e., the original TOMAS-30 default) is used 
for another set of simulations (denoted with “1 hr” added to 
the names of the simulations in Fig. 1). 

The TOMAS algorithm has two subroutines, condensa- 
tion/nucleation and coagulation, each of which uses an adap- 
tive internal time step. The 10 min or 1 h time step described 
here represents the “master” time step for microphysics that 
governs how frequently TOMAS alternates between conden- 
sation/nucleation and coagulation. In the case of the 1 h time 
step, condensation/nucleation occurs first for 1 h and then co- 
agulation occurs for the full hour. For the 10 min time step, 
six loops happen in the TOMAS microphysics algorithm 
(i.e., condensation/nucleation occurs first for 10 min, then co- 
agulation occurs for 10 min, back to condensation/nucleation 
for the next 10 min and so forth until the total time becomes 
1 h). Even when the “master” time step is 1 h, however, both 
condensation/nucleation and coagulation may subdivide this 
into smaller, internal time steps according to their adaptive 
schemes. 

The TOMAS-40 models with 10 min time steps are con- 
sidered to be the reference cases against which the other 
models are compared. Because of the uncertainties in the nu- 
cleation theories themselves and other challenges associated 
with global microphysical simulations, this does not imply 
that the reference cases are perfect simulations of the real 


Geosci. Model Dev., 6, 1221-1232, 2013 


www.geosci-model-dev.net/6/1221/2013/ 


Y. H. Lee et al.: Representation of nucleation mode microphysics 


1225 


(a) Binary— lnm Jnuc 




(d) Ternary— lnm Jnuc 
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Fig. 2. Pressure (hPa) vs. latitude maps for March-to-May averaged nucleation rates (/nuc) ( cm ^ s 'at 273 K and 1 atm), CN10 (Dp > 
10 nm) (cm - 3 at 273 K and 1 atm), CCN(0.2 %) (cm -3 at 273 K and 1 atm) for Binary-1 nm and Ternary-1 nm simulations. 


atmosphere. Rather, here we use the reference case to de- 
note a simulation with a numerically accurate simulation of 
nucleation mode dynamics and free of errors resulting from 
assumptions made in the Kerminen parameterization. 

To investigate the impact of uncertainties in nucleation 
rates on the representation of nucleation mode dynamics, 
each model configuration is used for two nucleation sce- 
narios: a low-nucleation-rate scenario with Vehkamaki et 
al. (2002) binary nucleation (nucleation is generally confined 
to the free troposphere) and a high-nucleation-rate scenario 
with Napari et al. (2002) ternary nucleation. Uncertainties in 
nucleation rates are represented using binary and ternary nu- 
cleation because the binary nucleation parameterization un- 
derpredicts nucleation rates in the boundary layer whereas 
the ternary nucleation parameterization produces unrealisti- 
cally high rates (Jung et al., 2006). We would like to note that 
the TOMAS model used in this study does not include any 
additional boundary nucleation scheme such as the cluster- 
activation nucleation scheme because it has been shown that 
the ternary nucleation scheme used here does predict fre- 
quent nucleation events in the boundary layer (Jung et al., 
2006, 2008, 2010). 

Figure 1 shows the ten simulations that will be discussed 
in this paper. The simulation names reflect the important fea- 
tures in each simulation: nucleation parameterization, size 
cutoff, and aerosol microphysical time step. As noted al- 
ready, the Binary- 1 nm and Ternary- 1 nm are the numerical 


benchmark simulations against which all biases in the aerosol 
number budgets are calculated. Regarding the computational 
burden of each model configuration, based on a one-month 
simulation with binary nucleation using a 600 Mhz single 
processor of an SGI Origin 300, it takes about 58 h for 
TOMAS-30. For TOMAS-36, it takes ~73 h using the 1 h 
time step and ~122h using the 10 min time step. Finally, 
for TOMAS-40, it takes ~ 79 h using the 1 h time step and 
~ 123 h using the 10 min time step. 


4 Results and discussion 

All model results presented in this paper are based on a three- 
month average from March to May using the GCM clima- 
tological meteorology after a three-month spin-up. Our re- 
sults are limited to three-month averages due to the com- 
putational demand in the TOMAS-36 and TOMAS-40 mod- 
els, but this time period should be sufficient for determining 
the relative errors due to differences in model-configuration 
assumptions. Figure 2 shows pressure-latitude maps of the 
three-month average nucleation rates (7 n uc)- CN10 (parti- 
cles with diameters larger than lOnm) concentrations, and 
CCN(0.2%) for the Binary- lnm and Ternary- lnm simu- 
lations. Figure 2 shows a very similar pattern to Fig. 3 in 
Pierce and Adams (2009c) that also presents annually aver- 
aged values of the same quantities. In Fig. 2a and d, the 7 nuc 
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Table 1. Global three-month average of the nucleation rate (./nuc), th e fonnation rate of particles larger than 10 nm from nucleated particles 
(J io), the ultrafine-mode (10 nm S Dp <10 nm) burden and lifetime, and the CCN mode (Dp > 70 nm) burden and lifetime in the six basic 
scenarios. Values nonnalized by tropospheric volume at 273 K and 1 atm, assuming 12 000 m as the height of the tropopause. 


Simulations 

•AlUC 

Tio 

Ultrafine mode 

Ultrafine mode 

CCN mode 

CCN mode 


[cm -3 day -1 ] 

[cm -3 day -1 ] 

burden [cm -3 ] 

lifetime [days] 

burden [cm -3 ] 

lifetime [days] 

Ternary- 1 nm 

3.4 x 10 8 

116 

758 

5.9 

142 

5.5 

Ternary-3 nm 

4.6 x 10 8 

126 

942 

7.0 

154 

5.7 

Ternary- 1 0 nm 

4.5 x 10 8 

3367 

1380 

0.4 

168 

5.7 

Binary- 1 nm 

4.2 x 10 4 

47 

458 

6.7 

115 

5.1 

Binary-3 nm 

1.2 x 10 4 

58 

538 

6.8 

120 

5.3 

Binary- 1 0 nm 

1.6 x 10 3 

365 

885 

2.4 

128 

5.5 


in Ternary- 1 nm is several orders of magnitude higher than 
that in Binary-lnm and is ubiquitous throughout the tropo- 
sphere in the regions that receive sunlight between March 
and May. The nucleation in Binary- 1 nm occurs in high alti- 
tudes and high latitude but does not predict much nucleation 
in the boundary layer. The distribution of CN10 concentra- 
tions, shown in Fig. 2b and e, is similar to the nucleation 
rate except for the surface layer where primary particles con- 
tribute significantly to CN10 concentrations. However, the 
difference in CN10 concentrations between the two nucle- 
ation mechanisms is much smaller than the difference in ./ nuc 
because a higher nucleation rate results in slower conden- 
sational growth rates (due to larger condensation sinks) and 
faster coagulation removal rates, which dampen the CN10 to 
changes in nucleation. Figure 2c and f show fairly insensi- 
tive CCN(0.2 %) to the nucleation rate because the probabil- 
ity of growth of nucleated particles to CCN-sized particles 
decreases with higher nucleation rate (Pierce and Adams, 
2009b). The similar results in Pierce and Adams (2009c) 
were annually averaged showing that these results are not 
largely sensitive to the choice of three-month seasons in the 
current work, and our key conclusions are unlikely to be 
changed by examining an annual average. 

4.1 Globally averaged aerosol number budgets 

Simulating aerosol microphysics in the nucleation mode 
with the Kerminen parameterization is influenced by the 
choice of the lowest size boundary. Table 1 presents the 
global-average aerosol budgets for the six simulation sce- 
narios including the J nuc (nucleation rate), the formation 
rate of particles larger than 10 nm from nucleated particles 
(7io), the number concentration and lifetime of ultrafine- 
mode particles (10nm< D p < 70 nm) and CCN mode par- 
ticles (70nm<Z)p): Binary-lnm, Ternary- 1 nm, Binary- 
3 nm, Ternary-3 nm, Binary-10 nm and Temary-10 nm. Com- 
pared to the lnm reference cases (i.e., Ternary-1 nm, Binary- 
1 nm), globally averaged 7 I111C is overpredicted using the 
3 nm and 10 nm boundaries for ternary nucleation by ~ 35 % 
but is underpredicted for binary nucleation by 70-96%. 
The Kerminen parameterization tends to overpredict globally 



Fig. 3. Pressure (hPa) vs. latitude maps for March-to-May averaged 
ratios of CN10 concentrations from the following scenarios: (a) 
Binary- lOnm to Binary- 1 nm, (b) Binary-3 nm to Binary- 1 nm, (c) 
Ternary- lOnm to Ternary- 1 nm, and (d) Ternary- 3 nm to Temary- 
1 nm. 


averaged Jio in the 3 nm and the 10 nm simulations. The 
overprediction of J io in the 3 nm simulations is only by 24 % 
for binary and 8% for ternary, whereas the lOnm simula- 
tions overpredict significantly by a factor of 7 for binary and 
a factor of 28 for ternary. Interestingly, with the Kerminen 
parameterization, the binary nucleation simulations overpre- 
dict iio despite underpredicting J auc . Compared to the refer- 
ence simulations, the aerosol burden and lifetime of the ul- 
trafine mode in the 3 nm simulations differ by a few percent 
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Table 2. Same as Table 1 but for the four sensitivity scenarios. 


Simulations 

^nuc 

4io 

Ultrafine mode 

Ultrafine mode 

CCN mode 

CCN mode 


[cm -3 day -1 ] 

[cm -3 day - 1 ] 

burden [cm -3 ] 

lifetime [days] 

burden [cm -3 ] 

lifetime [days] 

Ternary- 1 nm 

3.4 x 10 8 

116 

758 

5.9 

142 

5.5 

Ternary- 3 nm 

4.6 x 10 8 

126 

942 

7.0 

154 

5.7 

Binary- 1 nm 

4.2 x 10 4 

47 

458 

6.7 

115 

5.1 

Binary-3 nm 

1.2 x 10 4 

58 

538 

6.8 

120 

5.3 

Ternary- 1 nm-1 h 

3.3 x 10 8 

132 

753 

5.3 

148 

5.4 

Ternary- 3 nm-1 h 

4.4 x 10 8 

135 

942 

6.5 

160 

5.7 

Binary- 1 nm-1 h 

3.9 x 10 4 

50 

460 

6.4 

121 

5.2 

Binary-3 nm-1 h 

1.2 x 10 4 

62 

539 

6.5 

125 

5.3 


to ~ 25 %, and the 10 nm simulations show much larger er- 
rors than do the 3 nm simulations. Even though J \ o is over- 
predicted significantly in the 10 nm simulations, the aerosol 
burden in the ultrafine mode is overpredicted only by about 
80-90 %, which can be explained by the shorter lifetime of 
the ultrafine mode due to increased coagulation and reduced 
growth (keeping the particles at smaller sizes where coagu- 
lational losses are faster). In all cases, CCN budgets are less 
affected by the choice of lower size boundary than ultrafine- 
mode budgets (and thus CN10 budgets). 

We note that the underestimated 7 nuc in lOnm simulations 
(compared to the 1 nm simulations) for binary nucleation can 
be explained with the higher condensational sink (i.e., higher 
J\o). However, this is not the case for the ternary nucleation. 
During fast nucleation events by ternary nucleation, the small 
nucleation mode particles (1-3 nm) can be a relatively im- 
portant condensational sink as their number concentration is 
very high when nucleation occurs, resulting in lower nucle- 
ation rates. However, when nucleation is not occurring, the 
number concentration of the nucleation mode and the con- 
tribution of these particles to the condensational sink is gen- 
erally very small. This hypothesis can be also supported by 
the lower sulphuric acid concentration in TOMAS-40 over 
areas showing very fast nucleation rate. This explains why 
the J\f) overprediction is so much worse in the 10 nm sim- 
ulation (compared to the 1 nm simulation) with the ternary 
nucleation. 

In the GISS-TOMAS model, the Kerminen parameteriza- 
tion overpredicts 7io (or -h , defined as a particle formation 
rate for a particle larger than 3 nm) systematically, which 
leads to an overprediction of CN10 concentrations. Some 
of the known drawbacks in the Kerminen parameterization 
(e.g., no coagulational growth) do not explain the overpre- 
diction seen in the GISS-TOMAS model. Also, our results 
do not agree with Kerminen et al. (2004a), which does not 
show a consistent overprediction in J 3 , 7io, and total num- 
ber concentrations by the parameterization compared to their 
box-scale explicit dynamics model at same nucleation rates. 
The error in Jt, and J\o by the Kerminen parameterization 
might be improved with the revised formulation (Lehtinen 


et al., 2007; Anttila et al., 2010), but this is unlikely to ex- 
plain the overpredicted J\o seen in this work. Nevertheless, 
we found that the Kerminen parameterization with the 3 nm 
boundary introduces much less error in aerosol number pre- 
dictions than that with the 10 nm boundary, which is consis- 
tent with Kerminen et al. (2004a). 

Perhaps the most likely explanation for the overprediction 
of 7 10 seen here is that the Kerminen parameterization calcu- 
lates the survival probability from 1 mn to 3 or 10 nm based 
on the current conditions of the grid cell, and then the growth 
is assumed to occur instantaneously. This is the instantaneous 
growth rate assumption mentioned in Sect. 3. However, in re- 
ality, growth rates are frequently ~ 2 nm h - 1 or less so that 
growth to lOnm will take several hours during which con- 
ditions will change. Since photochemical activity, condens- 
able vapor concentrations, and growth rates are all quite high 
during the peak of a nucleation event, basing a calculated 
survival on instantaneous conditions during an event would 
likely lead to an overprediction in .7 1 0 . Moreover, the over- 
prediction would be expected to be much more severe with 
a 10 nm cutoff compared to a 3 nm cutoff, consistent with 
the results shown here. For a very pristine atmosphere with 
a nucleation event (i.e., convective cloud outflow areas), the 
overprediction in 7io may be worse by missing losses by self- 
coagulation within the nucleation mode by the Kerminen pa- 
rameterization. 

Unlike the lower size limit, global-average CN10 bud- 
gets are quite insensitive to the time step change in aerosol 
microphysical processes. Table 2 presents globally aver- 
aged aerosol number budgets for the remaining four simu- 
lations: Binary-1 nm-1 h, Ternary-1 nm-1 h, Binary-3 nm-1 h 
and Ternary-3 nm- 1 h. Compared to the simulations with a 
shorter time step presented in Table 1, aerosol budgets in 
these simulations show a small difference (mostly a few per- 
cent). Although small, we found some systematic differences 
when changing the aerosol microphysical time step. For ex- 
ample, with the 1 h time step, global-average 7 nuc is 2-6 % 
lower while global-average 7io is 6-13% higher. Global- 
average CCN mode burdens are higher by 4-5 %, which 
might be explained with 12 - 20 % higher condensational 
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3nm/ lnm 


1.08 


) Ternary lOnm/lnm 


1.64 


(d) Ternary 3nm/lnm 
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Fig. 4. Latitude vs. longitude maps of March-to-May average ratios 
of CN10 concentrations in the first vertical layer from the follow- 
ing scenarios: (a) Binary- 1 0 nm to Binary- 1 nm, (b) Binary-3 nm to 
Binary- lnm, (c) Ternary- lOnm to Ternary- 1 nm, and (d) Temary- 
3 nm to Ternary- 1 nm. The value on the top right of each figure is a 
global average of the ratios displayed in the map. 


growth rates for particles with diameters larger than 70 nm 
(not shown). Global-average CN10 burdens appear to be in- 
sensitive to the time step change. In conclusions, these results 
indicate that using a 1 h time step can lead to quite satisfac- 
tory results for these simulations but that a 10 nm lower size 
boundary can result in significant biases. 

4.2 Impact on spatial distributions of CN10 and 

CCN(0.2 %) 

Figure 3 shows pressure-latitude maps for March-to-May 
averaged ratios of CN10 concentrations in the 3 or lOnm 
simulations compared to those in the 1 nm benchmark sim- 
ulations. Unlike global-average quantities, they show much 
larger differences in CN10 in some regions with the maxi- 
mum being more than a factor of four. The overpredictions 
of CN10 in the 3 nm and 10 mn simulations occur throughout 
the troposphere for both binary and ternary nucleation, and 
the large differences occur at high-nucleation-rate regions as 
shown in Fig. 2. The higher CN10 concentrations in the 3 
and 10 nm simulations are likely a result of the overpredicted 
J io or by the Kerminen parameterization, although other 

factors may play a role. Unlike the 3 nm simulations, the 
10 nm simulations (shown in Fig. 3a and c) show underpre- 
dicted CN10 in some parts of the tropics for both binary and 
ternary nucleation schemes. The underpredicted CN10 areas 
show no nucleation events for the binary nucleation and, for 
the ternary nucleation, there are noticeably lower nucleation 
rates and lower free-ammonia concentrations in TOMAS-30 
than TOMAS-40. This clearly shows that the choice of the 
lowest size boundary is important for CN10 prediction. Sim- 
ilar to the global-average CN10 predictions, the zonal aver- 
age of ratios of CN10 concentrations in the 3 nm simulations 


(b) Binary 3nm/lnm 




(d) Ternary 3nm/lnm 



-90 -60 -30 0 30 60 90 


.9 .95 .99 1.01 1.05 1.1 


Fig. 5. Pressure (hPa) vs. latitude maps for March-to-May averaged 
ratios of CCN(0.2%) concentrations from the following scenar- 
ios: (a) Binary-10 nm to Binary-1 nm, (b) Binary-3 nm to Binary- 
1 nm, (c) Ternary- lOnm to Ternary- 1 nm, and (d) Ternary-3 nm to 
Ternary- 1 nm. 


to those in the 1 nm simulations with 1 h time steps are al- 
most identical to the 10 min time step shown in Fig. 3b and 
d. The change in the CN10 predictions when the time step 
is changed from the 10 min to 1 h is only a few percent (not 
shown). 

Figure 4 shows latitude-longitude maps of ratios of the 
lowermost-layer CN10 concentrations in the 10 nm or 3 nm 
simulations compared to those in the 1 nm benchmark simu- 
lations. The spatial patterns of the CN10 prediction errors are 
different between the two nucleation mechanisms. Compar- 
ing lower size cutoffs, the patterns are similar but the lOnm 
cutoff consistently shows larger errors. Averaged across the 
entire lowermost layer, the 3nm simulations introduces ~ 10 
to ~ 20 % difference for the binary and ternary, respectively; 
in the 10 nm simulations, the differences are ~ 15 to ~ 60 %. 
The surface-layer CN10 in the 3 nm simulations agrees with 
the 1 nm benchmark simulation within 30 % in most regions, 
but the error can be a factor of 2 in the Middle East, where 
there is a high nucleation rate in Ternary-3 nm, reflecting the 
overprediction of J\o in the 3 nm simulation. Again, for the 
3 nm and 1 nm simulations with the 1 h time step, the surface- 
layer CN10 agrees very well (within a few percent in most 
regions) to those with 10 min time step (not shown). The 
small impact on CN10 by the choice of time step may be a 
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Fig. 6. Latitude vs. longitude maps of March-to-May average ratios 
of CCN(0.2 %) concentrations in the first vertical layer from the 
following scenarios: (a) Binary-10 nm to Binary-1 nm, (b) Binary- 
3nm to Binary- lnm, (c) Ternary- lOnm to Ternary- lnm, and 
(d) Ternary-3 nm to Ternary- 1 nm. The value on the top right of 
each figure is a global average of the ratios displayed in the map. 


result of higher condensational growth rate and shorter life- 
time (Sect. 4.1). 

Figure 5 shows zonal-average ratios of CCN(0.2 %) in the 
3 nm and 10 nm simulations with those in the 1 nm bench- 
mark simulations. Unlike Fig. 3, the spatial distribution of 
CCN(0.2 %) ratios is similar between binary and ternary nu- 
cleation, but the 3nm and lOnm simulations are somewhat 
different from each other. Overall error in the CCN(0.2 %) 
prediction by using a 3 nm cutoff is no more than about 10 % 
in any location. For the 10 nm cutoff, the error is generally 
less than 20%, but is as high as 50% in some locations. 

The spatial distribution of the lowermost-layer 
CCN(0.2%) ratios (shown in Fig. 6) is similar to, but 
smaller than, the lowermost-layer CN10 (shown in Fig. 4). 
The CCN(0.2%) in the 3 nm simulations differ from the 
1 nm benchmark by less than 5 % in most regions, while 
that in the 10 nm simulations is overpredicted by 10-50% 
in most of the Northern Flemisphere for both nucleation 
mechanisms. Overall, CCN(0.2 %) is affected by the choice 
of lower size limit but to a much lesser degree than CN10, 
and the CCN(0.2%) difference is quite similar between 
two nucleation schemes, which reflects the insensitivity of 
CCN(0.2 %) to nucleation. 

The choice of microphysical time step has a negligible 
impact on CN10 but a greater influence on CCN(0.2%), at 
least in some locations. Figure 7 presents spatial distribu- 
tions of CCN(0.2 %) ratios of Ternary-3 nm-1 h to Ternary - 
3 nm and Ternary-1 nm-1 h to Ternary-1 nm simulations. The 
binary nucleation cases are not shown because they have 
similar spatial trends and magnitudes. The biases in the 
CCN(0.2 %) using the 1 h time step in both zonal averages, 
shown in Fig. 7, are quite similar between the 3 nm and the 
1 nm cutoffs and are mostly within 5-10%; the same is true 



Fig. 7. Same as Fig. 5 but for the following scenarios: (a) Temary- 
3 nm- 1 h to Temary-3 nm and (b) Ternary- 1 nm- 1 h to Ternary- 1 nm. 


at the model surface (not shown). In fact, when comparing 
with the 1 nm simulation with 1 h time step, the biases in 
the CCN(0.2%) from the lOnm simulations turn out to be 
smaller (see Figs. 8 and 9) than the comparisons with the 
benchmark simulations (see Figs. 5 and 6) but still larger 
than the 3 mn simulations with the 1 h time step. Also, Fig. 8 
shows a very similar spatial pattern of the CCN(0.2 %) errors 
between the 10 nm and 3 nm simulations. Averaged across 
the entire lowermost layer, using the same time step, the 
lOnm simulations introduces a ~3 to ~6% difference for 
the binary and ternary, respectively (see Fig. 9); using dif- 
ferent time steps, the errors by the lOnm simulations are 
7% and 9% (see Fig. 6). This indicates that the larger 
CCN(0.2%) deviation by the 10 nm simulations (shown in 
Fig. 5a and c) is partly attributed to the different time step. 
Thus, the choice of a time step is as important in predicting 
CCN(0.2 %) as the choice of the lower size limit. 


5 Conclusions 

We investigated the representation of the nucleation mode 
(1 nm < D p < 10 nm) particle dynamics in a global model 
with aerosol microphysics by comparing explicit represen- 
tations of nucleation mode aerosol to parameterizations of 
nucleation mode microphysics. This study uses the global 
aerosol microphysics model GISS-TOMAS, varying its low- 
est aerosol diameter boundary: 1 nm, 3 nm, and 10 nm. The 
microphysics of nucleation mode particles are explicitly re- 
solved with the 1 nm boundary. The model with the lOnm 
and 3 nm boundaries uses a nucleation mode dynamics pa- 
rameterization proposed by Kerminen et al. (2004a), to ac- 
count for the growth of nucleated particles to lOnm and 
3 nm, respectively. We also compared a 10 min time step ver- 
sus a 1 h time step in the TOMAS algorithm to investigate 
the impact of time step on aerosol number predictions. The 
simulations with 1 nm size cutoff and 10 min time step are 
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Fig. 8. Same as Fig. 5 but for the following scenarios: (a) Binary- 
10 nm to Binary- 1 nm-1 h, (b) Binary-3 nm-1 h to Binary- 1 nm-1 h, 
(c) Ternary- lOnm to Ternary- 1 nm-1 h, and (d) Ternary-3 nm-1 h to 
Ternary- 1 nm-1 h. 


used as reference cases to estimate errors caused by increas- 
ing lower size cutoff and time step. 

Different representations of the nucleation mode have a 
minor effect on globally averaged CCN mode burdens and 
lifetimes. However, they do affect global-average 7io (forma- 
tion rate of particles greater than lOnm from nucleated par- 
ticles) and the lifetime and burden of ultrafine particles (and 
thus CN10). When the lower cutoff diameter is raised from 
1 nm to 3 nm and 10 nm, it leads to systematic biases in 7 io, 
CN10 and CCN(0.2 %), and these biases are generally much 
greater with the lOnm cutoff. The CN10 concentrations are 
biased high, which may be caused by overpredictions in J 3 
or 7io by the Kerminen parameterization. The overpredicted 
7io and J 3 are likely the result of the instantaneous growth 
rate assumption rather than missing coagulational growth. 

The magnitude of the errors in predicted CN10 concentra- 
tions depends strongly on location and what nucleation pa- 
rameterization is used. Not surprisingly, simple nucleation 
microphysics leads to larger errors in regions with stronger 
nucleation and when using the faster ternary nucleation rates. 
The lOnm and 3 nm boundary simulations show errors in 
zonally averaged CN10 predictions up to a factor of 3-5 
where high nucleation rates occur. However, the lowermost- 
layer CN10 deviations are mostly within 50% except for 
ternary nucleation with 10 nm simulations. The change in the 


(a) Binary lOnm/lnm- lhr 1.03 (b) Binary 3nm-lhr/lnm- lhr 1.01 
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Fig. 9. Same as Fig. 6 but for the following scenarios: (a) Binary- 
10 nm to Binary- 1 nm-1 h, (b) Binary-3 nm-1 h to Binary- 1 nm-1 h, 
(c) Ternary- lOnm to Ternary- 1 nm- 1 h, and (d) Ternary- 3 nm-1 h 
to Ternary- 1 nm- 1 h. The value on the top right of each figure is 
a global average of the ratios displayed in the map. 


aerosol microphysics time step from 10 min to 1 h has little 
influence on CN10 budgets, possibly because of the compen- 
sation between the higher 7io and the short lifetimes. 

The errors in zonal-average CCN(0.2 %) with the 3 nm 
boundary are mostly within 10% of the 1 nm boundary 
case, while those for the lOnm boundary case are larger 
(within 20% in most regions). Similarly, the surface-layer 
CCN(0.2 %) from the 3 nm boundary cases agree with the 
benchmark 1 nm model as most regions have only 1 % to 5 % 
differences, but the CCN(0.2%) from the 10 nm boundary 
cases differ by more than 10-20% in most northern hemi- 
spheric areas. We found that this larger deviation is partly 
attributed to the time step, which is different between the 
benchmark simulations (i.e., 10 min) and the 10 nm simula- 
tions (i.e., 1 h). Comparing the simulations with the 1 h time 
step, the deviations in the CCN(0.2%) prediction with the 
10 nm boundary are reduced noticeably but still larger than 
those with the 3 nm boundary. 

This study shows that the representation of nucleation 
mode in a global aerosol microphysics model has a minor in- 
fluence on CCN(0.2 %) overall, and by extension the aerosol 
indirect effect, but a more significant impact on CN10 con- 
centrations. For CCN simulations, a model with a lower di- 
ameter limit of 1 0 mn can be sufficient, but the 3 nm limit is 
recommended. Studies focused on CN10 predictions or nu- 
cleation will benefit from using a 1 nm lower size boundary. 
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